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Abstract 



We present a proposal for performing model-independent jets plus missing energy 
searches. Currently, these searches are optimized for mSUGRA and are consequently 
not sensitive to all kinematically-accessible regions of parameter space. We show that 
the reach of these searches can be broadened by setting limits on the differential 
cross section as a function of the total visible energy and the missing energy. These 
measurements only require knowledge of the relevant Standard Model backgrounds 
and can be subsequently used to limit any theoretical model of new physics. We 
apply this approach to an example where gluinos are pair-produced and decay to 
the LSP through a single-step cascade, and show how sensitivity to different gluino 
masses is altered by the presence of the decay chain. The analysis is closely based 
upon the current searches done at the Tevatron and our proposal requires only small 
modifications to the existing techniques. We find that within the MSSM, the gluino 
can be as light as 125 GeV. The same techniques are applicable to jets and missing 
energy searches at the Large Hadron Collider. 



I. INTRODUCTION 



One of the most promising signatures for new physics at hadron coUiders are events with 
jets and large missing transverse energy (^r)- These searches are very general and cover 
a wide breadth of potential new theories beyond the Standard Model. Jets + searches 
pose a significant challenge, however, because the Standard Model background is difficult 
to calculate in this purely hadronic state. The general nature of the signature motivates 
performing a search that only requires calculating the Standard Model background. The 
challenge, then, is to minimize the risk of missing new physics while still accounting for 
our limited understanding of the background. All experimental searches of jets + at 
hadron colliders have been model-dependent, attempting to be sensitive to specific models 
[H El [31 m El E] . Initial studies for the Large Hadron Collider (LHC) have been dominantly 
model-dependent Ei j^. In this article, we explore how modest modifications to the 
existing jets and studies can allow them to be model-independent, broadening the reach 
of the experimental results in constraining theoretical models. 

Currently, jets plus searches at the Tevatron are based on the Minimal Supersym- 
metric Standard Model (MSSM) ^Hj and look for production of gluinos {g) and squarks 
(g), the supersymmetric partners of gluons and quarks, respectively [21 El [4]. These par- 
ticles subsequently decay into the stable, lightest supersymmetric particle (LSP), which is 
frequently the bino, the supersymmetric partner of the photon. The MSSM contains hun- 
dreds of parameters and it is challenging to place mass bounds in such a multi-parameter 
space. To make this tractable, the CMSSM (or mSUGRA) ansatz has been used fT2]. The 
CMSSM requires common scalar masses (mo), gaugino masses (mi ), and trilinear scalar soft 
couplings (^40) at the unification scale, in addition to electroweak symmetry breaking, gauge 
coupling unification, and R-parity conservation. The entire particle spectrum is determined 
by five parameters. 

One important consequence of this theory is that the ratio of gaugino masses is fixed 
at approximately rrig : m^ : m^ 6 : 2 : 1, where W refers to the triplet of winos 

{W^^ W^)^ the supersymmetric partners of the electroweak gauge bosons. Due to the number 
of constraints in the CMSSM, the bino is the LSP throughout the range of parameter space 
that the Tevatron has access to. Furthermore, due to the renormalization group running 
of the squark masses, the squarks are never significantly lighter than the gluino. Thus, 
the ratio in masses between the lightest colored particle and the LSP is essentially fixed. 
The CMSSM is certainly not representative of all supersymmetric models (see, for example, 
p!3l [T4[ ITSl [T6[ [TTl [T8]). let alone the wider class of beyond the Standard Model theories 
that jets and searches should have sensitivity to. Verifying that a jets and search 
has sensitivity to the CMSSM does not mean that the search is sensitive to a more generic 
MSSM. 

Existing searches for gluinos and squarks make strong assumptions about the spectrum 
and it is unclear what the existing limits on squark-like and gluino-like particles are. Because 
squarks have electric charge, LEP can place limits of 92 GeV on their mass [19j; however, 
gluinos do not couple to either the photon or and so limits from LEP2 are not strong. Cur- 
rently, the tightest model-independent bound on color octet fermions (such as gluinos) comes 
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from thrust data at ALEPH [20] and OPAL [2T]. New colored particles should contribute 
at loop-level to the running of the strong coupling constant ag. To date, the theoretical 
uncertainties in the value of ag have decreased its sensitivity to new particle thresholds. Ad- 
vances in Soft-Collinear Effective Theory, however, have been used to significantly reduce the 
uncertainties in ag from LEP data. The current bound on color octet fermions is 51.0 GeV 
at 95% confidence [22j; no limit can be set for scalar color octets. 

There is no unique leading candidate for physics beyond the Standard Model; therefore, 
searches for new physics need to be performed in many diflFerent channels. Ideally, one 
should perform totally model-independent searches that only employ the Standard Model 
production cross section for physics with the desired channels and the correct kinematics. 
The goal is to be sensitive to a large number of different models at the same time so that 
effort is not wasted in excluding the same parts of Standard Model phase space multiple 
times. 

Some progress on experimental model-independent searches has been made. In an ambi- 
tious program, the CDF Collaboration at the Tevatron has looked at all possible new channels 
simultaneously (i.e.. Vista, Sleuth, Bumphunter) [23", "24^, "25 ] ; however, these searches have 
some drawbacks over more traditional, channel-specific searches. The most important draw- 
back is that it is difficult, in the absence of a discovery, to determine what parts of a given 
model's parameter space are excluded. 

On the theoretical front, MARMOSET [26] is a hybrid philosophy that attempts to 
bridge model-independent and model-dependent searches with the use of On-Shell Effective 
Theories (OSETs). OSETs parameterize the most experimentally relevant details of a given 
model - i.e., the particle content, the masses of the particles, and the branching ratios of 
the decays. By using an on-shell effective theory, it is possible to easily search through all 
experimentally relevant parameters quickly. The on-shell approximation is not applicable in 
all situations, but OSETs can still give a rough idea of where new physics lies. 

In this article, we will explore the discovery potential of jets and missing energy chan- 
nels. In previous work [28], we presented a simple effective field theory that can be used to 
set limits on the most relevant parameters for jets and missing energy searches: the masses 
of the particles. While this approach seems obvious, existing searches at hadron colliders 
(Tevatron Run II, Tevatron Run I, UA2, UAl) are based on CMSSM-parameterized super- 
symmetry breaking. The previous paper studied how varying the decay kinematics changed 
the sensitivity of the searches and pointed out regions of parameter space where sensitivity 
is particularly low due to kinematics. However, this gluino-bino module was still a model- 
dependent analysis in that it assumed pair-production of a new colored fermionic particle 
directly decaying to a fermionic LSP. 

This paper will extend the analysis in two ways. First, we propose a completely model- 
independent analysis for jets and missing energy searches. This approach only requires 
knowledge of the Standard Model and places limits on differential cross sections, from which 
it is possible to set model-dependent limits. In the second portion of the paper, we use this 
approach to extend our previous analysis of a directly decaying colored particle to contain 
a single-step cascade and study how this altered spectrum affects the final limits on the 
gluino's mass. 
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II. OVERVIEW OF MODELS 



Before continuing with the main theme of the article, let us take a moment to describe 
the class of models that jets + searches are sensitive to. There are two general classes 
of particle spectra that will be covered by such searches, each of which has a stable neutral 
particle at the bottom of the spectrum. Typically, the stability of these neutral particles is 
protected by a discrete symmetry (e.g., R-parity, T-parity, or KK-parity) and, consequently, 
these particles are good candidates for the dark matter. In one class of models, the theory 
contains a new colored particle that cascade decays into the dark matter. In the other 
class, new elect roweak gauge bosons are produced. The dark matter particle may either be 
produced along with the new bosons, or may be the final step in their decays. 

The first class can be thought of as being generally SUSY-like where the lightest colored 
particle is dominantly produced through the Standard Model's strong force. The lightest 
colored particle then cascade decays down to the stable, neutral particle at the bottom of 
that sector. These cascades will either be lepton-poor or lepton-rich. Lepton-poor cascades 
occur when there is no state accessible in the cascades that have explicit lepton number (e.g., 
sleptons) and frequently occur when the cascades are mediated by W^^ Z^, or Higgs bosons. 
A simple supersymmetric example of a lepton-poor cascade decay is a theory where the scalar 
masses are made heavy and only gauginos and Higgsinos are available in the decay chains. 
This occurs, for instance, in PeV supersymmetry models, where the scalars are around 1000 
TeV and the fermions of the MSSM are in the 100 GeV to 1 TeV range. Producing the color- 
neutral states of such a theory is difficult at hadron machines; consequently, the production 
of new particles will occur primarily through the decay of the gluino. 

One potential cascade decay of the gluino, which will be considered in further detail in 
the second half of the paper, is 

g qiq2W qiq2qzqAB. (1) 

In this cascade, the W decays directly into the B and a W^^ boson, which subsequently 
decays to two jets. This single-step decay is the dominant cascade if the gaugino masses 
are unified at high energies; in this case, the branching ratio of the gluino into the wino 
is ^ 80%. While these cascade decays are to some degree representative, the precise mass 
ratio of rUg : : makes a significant difference in the searches. In the limit where 
—> the energy from q^ and q^ is small, while if rUg the jets from qi and q2 

are soft. If > m^, this cascade is forbidden. Interestingly, spectra with unified gaugino 
masses are the most difficult to see because all four jets are fairly hard and diminish the 
missing energy in the event in comparison to the direct decay of the gluino, g qiq2B. 
Leptons from the decay of the W^^ boson can be used in the analysis as well (see 



Sec. YE). However, jets + + lepton studies are better suited for lepton-rich cascades. 
The addition of leptons to the searches makes the experimental systematics easier to control 
and improves trigger efficiencies. Not all spectra of new physics can be probed with these 
types of searches, though, and they are thus complimentary to the jets + search. 

Other cascades may produce a greater number of jets as compared to ([T]). In NMSSM 
theories where there is a new singlino at the bottom of the spectrum [27j, it is possible to 
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have cascade decays that start with the gluino, go to wino plus two jets, then bino plus two 
additional jets, and conclude with the singlino plus two more jets. The additional step in the 
decay process further diminishes the amount of missing energy in typical events, resulting in 
reduced limits on spectra. Other models, such as Universal Extra Dimensions (UEDs) 
and Little Higgs models with T-parity [30] also have new colored particles that subsequently 
cascade decay. The details of the exact spectra can alter the signal significantly as jets can 
become soft and missing energy is turned into visible energy. 

It is also possible that new electroweak gauge bosons are produced, which then cascade 
decay, producing jets before ending with the neutral stable particle. Little Higgs models 
with T-parity are one such example. In such models, the new heavy bosons and are 
produced through s-channel processes. The can decay to the and the dark matter 
Ah^ while the can decay to the Ah and higgs. It is also possible to produce the 
directly with the Ah through an s-channel boson. This vertex, however, is suppressed 
in comparison to the other two. 

III. PROPOSED ANALYSIS STRATEGY 

At the Tevatron, the jets + channel is divided into four separate searches (monojet, 
dijet, threejet, and multijet), with each search defined by jet cuts (9(30 GeV). Cuts on the 
missing transverse energy and total visible energy^ Ht of each event take place during the 
final round of selection cuts. The and Ht cuts are optimized for "representative" points 
in CMSSM parameter space for each of the (inclusive) Ij — 4:'^j searches. However, these 
and Ht cuts may not be appropriate for theories other than the CMSSM. Indeed, considering 
the full range of kinematically allowed phase space means accounting for many combinations 
of missing and visible energy. A set of static cuts on ^t and Ht is overly-restrictive and 
excludes regions of phase space that are kinematically allowed. 

This is explicitly illustrated in Fig. [T} which shows the ^t distribution of a dijet sample 
passed through two diflFerent sets of ^t and Ht cuts. The signal, a 210 GeV gluino directly 
decaying (i.e., no cascade) to a 100 GeV bino, is shown in white and the Standard Model 
background, in gray. The plot on the left shows the events that survive a 300 GeV Ht cut. 
While the Ht cut significantly reduces the background, it also destroys the signal above 
the ^T cut of 225 GeV. These cuts were used in the dijet search; they are optimized 
for a ^ 400 GeV gluino, but are clearly not ideal for the signal point shown here. A more 
optimal choice of cuts is shown on the right. While the lower Ht cut of 150 GeV keeps more 
background, it also keeps enough signal for a reasonable S/B ratio at low ^t- Therefore, 
with a ^T cut of 100 GeV, exclusion limits on this point in parameter space can be placed. 

A model-independent search should have broad acceptances over a wide range of kine- 
matical parameter space. Ideally, searches should be sensitive to all possible kinematics by 
considering all appropriate ^t and Ht cuts. This can be eflFectively done by plotting the 



The total visible energy Ht is defined as the scalar sum of the transverse momenta of each jet. 
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FIG. 1: Comparison of cuts and optimized cuts for a sample dijet signal for rrig = 210 GeV 
and = 100 GeV. Background distribution is shown in gray and signal distribution in white. 
(Left) Using the cuts Ht > 300 GeV and > 225 GeV (Right) Using the more optimal cuts 
Ht > 150 GeV and ^ 100 GeV. The optimized cuts allow us to probe regions with larger S/B. 



differential cross section as a function of and Ht^ 

(fa 



dHTd^T 



AHtA^t- (2) 



In this case, the results of a search would be summarized in a grid, where each box contains 
the measured cross section within a particular interval of and Ht- 

As an example, the differential cross section grids for exclusive Ij — 4+^' searches (see 
Table [l] for jet selection criteria) at the Tevatron are shown in Table IL The grids are made 
for the Standard Model background, which include W^ + nj^ Z^ + nj^ and tt + nj. The QCD 
background was not simulated; we expect the QCD contributions to be important for points 
in the lowest bin. For details concerning the Monte Carlo generation of the backgrounds, 
see Sec. HVBl 

From these results, it is straightforward to obtain limits on the differential cross section 
for any new physics signal. Consider a specific differential cross section measurement that 
measures A^^ events in an experiment. The Standard Model predicts B events, while some 
specific theory predicts B + S events, where S is the number of signal events. 

The probability of measuring n events is given by the Poisson distribution with mean 
fi = B + S. The mean /x is excluded to 84% such that 



nl 



excl (/X^^^l) 



E^<«'-16^ (3) 



n=0 

The solution to this equation gives the excluded number of signal events 

S'^^''\N^,B)=fi'-''\Nj-B. (4) 
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TABLE I: Summary of the selection criteria for the four exclusive (i.e., non-overlapping) searches. 
The two hardest jets are required to be central (I77I < 0.8). All other jets must have I77I < 2.5. 



The expected limit on the signal is then given by 



Nrr 



Nm=0 



(5) 



In the limit of large 5, the probability distribution approaches a Gaussian and we expect 
that 

lim (5^"^^ (5)) = (6) 

B^oo 



In the limit of small 5, we expect that 



\iml,S^''''\B)) = -ln(0.16) ^ 1.8. 



(7) 



The right column of Table [TT] shows the limit on the differential cross section for any 
new physics process. When presented in this fashion, the experimental limits are model- 
independent and versatile. With these limits on the differential cross section, anyone can 
compute the cross section for a specific model and make exclusion plots using just the signal 
limits shown in Table [TH For the comparison to be reliable, the detector simulator should be 
properly calibrated. 

In addition to the statistical uncertainty, systematic uncertainties can also be important. 
Unlike the statistical uncertainties, the systematic uncertainties can be correlated with each 
other. One important theoretical uncertainty is the higher-order QCD correction to the 
backgrounds. These QCD uncertainties result in K-factors that change the normalization 
of the background, but do not significantly alter the background shapes with respect to Ht 
and Because this uncertainty is highly correlated between different differential cross 
section measurements, treating the uncertainty as uncorrelated reduces the sensitivity of the 
searches. If a signal changes the shape of the differential cross section, e.g. causing a peak in 
the distribution, higher order corrections would be unlikely to explain it. To make full use of 
the independent differential cross section measurements, a complete error correlation matrix 
should be used. In practice, because the backgrounds are steeply falling with respect to 
Ht and ^t-> assigning an uncorrelated systematic uncertainty does not significantly hurt the 
resolving power of the experiment. In Table [TTl we have assigned a systematic uncertainty 
of egys = 50% to each measurement, which should be added in quadrature to the statistical 
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TABLE IL Differential cross section (in fb) for the Standard Model background is shown in the left 
column for exclusive Ij — 4+ j searches. The expected signal sensitivity at 84% confidence is shown 
on the right (in fb). The statistical error is shown to the left of the © and the systematic error is on 
the right. For purposes of illustration, we assume a 50% systematic error on the background. The 
gray boxes are kinematically forbidden. These results are for 4 fb~^ luminosity at the Tevatron. 



uncertainty. This roughly corresponds to the requirement that the total signal to background 
ratio is one. 
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The reduced chi-squared Xn value for measurements is 



where is the number of signal events and Bj is the number of background events in the 
j^^ box of the grid. The statistical error SLj and the systematic error Cgys x Bj is read off 
from Table [11} In order to have a useful significance limit, it is necessary to only include 
measurements where there is an expectation of statistical significance; otherwise, the x% 
diluted by a large number of irrelevant measurements. There is no canonical way of dealing 
with this elementary statistical question, although the CL^- method is the most commonly 
used [311 132]. In this article, we take a very simple approach. If the expected significance 
for a single measurement is greater than a critical number, S'^^^*, it is included in the Xn^ 
otherwise it is not. We tried several values of S^^^^ and the experimental sensitivity to different 
theories was not altered by the different choices. We chose S^^^^ = 0.5 for the exclusion plots. 
This method does not maximize the reach in all cases, but because there are usually just 
a few measurements that give large significance, we are relatively insensitive to the exact 
statistical procedure. 

In what follows, we will apply the general philosophy presented here to find the exclusion 



region for gluinos that are pair-produced at the Tevatron.^ In Sec. IV, we will explain how 
the signal and background events have been generated. In Sec. [V| we will show how mass 
bounds can be placed on the gluino and bino masses using the proposed model-independent 
analysis and will discuss the challenges presented by cascade decays. We conclude in Sec. 

lYn 



IV. EVENT GENERATION 



A. Signal 

In this section, we discuss the generation of signal events for the gluino cascade decay 
shown in ([T]). The experimental signatures of this decay chain are determined primarily 
by the spectrum of particle masses. In particular, the mass splittings determine how much 
energy goes into the jets as opposed to the bino - i.e., the ratio of the visible energy to 
missing transverse energy. Events with large Ht and will be the easiest to detect; this 
is expected, for example, when a heavy gluino decays into a wino that is nearly degenerate 
with either the gluino or the bino. The reach of the searches is degraded, however, when the 
wino is included as an intermediate state in the decay chain. When the jets from the cascade 
decay are all hard, the missing energy is significantly smaller than what it would be for the 
direct decay case. Picking out signals with small missing transverse energy is challenging 



^ Throughout this article, "gluino" refers to a color octet fermion, "wino" to a charged SU(2) fermion, and 
"bino" to a neutral singlet. These names imply nothing more than a particle's quantum numbers. 
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because they push us closer to a region where the dominant background is coming from QCD 
and is poorly understood. This happens, in particular, when the mass splitting between the 
gluino and bino is large and the wino mass is sufficiently separated from both. When the 
wino is nearly degenerate with either the gluino or the bino, then we expect to see 2 hard jets 
and 2 soft jets from the decay. This case begins to resemble the direct decay scenario; there 
is more missing energy and, therefore, the signal is easier to see. It is particularly challenging 
to probe regions of parameter space where the gluino is nearly degenerate with the bino. For 
this case, even in the light-gluino region {rrig < 200 GeV), the benefit of the high production 
cross section for the gluinos is overwhelmed by the small missing transverse momentum in 
each event; the jets in these events are soft and the pt of the two binos approximately cancel 
when summed together [28j. Even if the gluinos are produced at large invariant mass, the 
situation is not markedly improved; in this case, the jets from each gluino are collinear and 
aligned with the ^t- Such events are easily mistaken as QCD events and eliminated by the 
cuts that are implemented to reduce the QCD backgrounds. 

The inclusion of hard initial-state jets significantly increases the exclusion reach in this 
degenerate region of parameter space. The initial-state radiation boosts the gluinos in the 
same direction, decreasing the angle between them, which in turn, enhances the ^t- There- 
fore, ISR jets allow us to capitalize on the high production cross section of light gluinos to 
set bounds on their masses. 

To properly account for initial-state radiation (ISR) and final-state radiation (FSR), 
MadCraph/MadEvent [33j was used to generate events of the form 

pp^ gg + Nj, (9) 

where = 0,1,2 is the multiplicity of QCD jets. Pythia 6.4 [34j was used for parton 
showering and hadronization. Properly counting the number of events after parton showering 
requires some care. In general, an (n + l)-jet event can be obtained in two ways: by a (n + 1) 
hard matrix-element, or by hard radiation emitted from an n-parton event during showering. 
It is important to understand which of the two mechanisms generates the (n + l)-jet final 
state to ensure that events are not double-counted. 

In this article, a version of the so-called MLM matching procedure implemented in Mad- 
Craph/MadEvent and Pythia |35j was used for properly merging the different parton multi- 
plicity samples. This matching has been implemented both for Standard Model production 
and for beyond the Standard Model processes. In this procedure, parton-level events are 
generated with a matrix element generator with a minimum distance between partons char- 
acterized by the k± jet measure: 

= Ai?|min(p|.,p|^.) 
beam) = p^-^ (10) 

where Ai?^^- 2[cosh(A7y) — cos(A(/))] [36]. The event is clustered using the clustering 
algorithm, allowing only for clusterings consistent with diagrams in the matrix element, 
which can be done since MadCraph generates all diagrams for the process. The (P values 
for the different clustered vertices are then used as scales in the as value corresponding to 
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that vertex, i.e. the event weight is multiphed by Yii ) ^ ^^^^^ the product is over the 
clustered vertices i. This is done in order to treat radiation modeled by the matrix element 
as similarly as possible to that modeled by the part on shower, as well as to correctly include 
a tower of next-to-leading log terms. A minimum cutoff d{i^j) > Qmin placed on all the 
matrix-element multi-part on events. 

After showering, the partons are clustered into jets using the standard k± algorithm. 
Then, the jet closest to the hardest parton in (ry, (/))-space is selected. If the separation 
between the jet and parton is within some maximum distance, rf(parton, jet) < Q^^^^ the jet 
is considered matched. The process is repeated for all other jets in the event. In this way, 
each jet is matched to the parton it originated from before showering. If an event contains 
unmatched jets, it is discarded, unless it is the highest multiplicity sample. In this case, 
events with additional jets are kept, provided the additional jets are softer than the softest 
parton, since there is no higher-multiplicity matrix element that can produce such events. 
The matching procedure ensures that jets are not double-counted between different parton 
multiplicity matrix elements, and should furthermore give smooth differential distributions 
for all jet observables. The results should not be sensitive to the particular values of the 
matching parameters, as long as they are chosen in a region where the parton shower is a 
valid description. Typically, the matching parameters should be on the order of the jet cuts 
employed and be far below the factorization scale of the process. For the gluino production, 
the parameters were 

Q^^ = 20 GeV = 30 GeV. (11) 

c 

Q-10 



1 

10-^ 

10-2 



Differential Jet Rate 1 

FIG. 2: Differential — > 1 jet rate for a matched sample of light gluino production. The full black 
curve shows the matched distribution, and the broken curves show the contributions from different 
matrix element parton multiplicity samples. The matching scale Q^f^ is marked by the dashed 
line. The full red curve shows the result using Pythia only. 
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Figure [2] shows the differential jet rate going from zero to one jets D{lj Oj), which is 
the maximum k± distance for which a Ij event is characterized as a Oj event. Below Qmfn? 
all jets come from part on showering of the Oj multiplicity sample. Above Qm^n^ the jets 
come from initial-state radiation. The main contributions in this region are from the Ij and 
2j multiplicity samples. The sum of all the multiplicity samples is a smooth distribution, 
eliminating double counting between the different samples. 

The simulations were done using the CTEQ6L1 PDF and with the renormalization and 
factorization scales set to the gluino mass [37j . The matched cross-sections were rescaled to 
the next-to-leading-order (NLO) cross sections obtained using Prospino 2.0. PGS was used 
for detector simulation [39], with jets being clustered according to the cone algorithm, with 
AR = 0.5. As a check on this procedure, we compared our results to the signal point given 
in [3j and found that they agreed to within 10%. 
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FIG. 3: Comparison of matched and unmatched events for a dijet sample of 150 GeV gluinos 
directly decaying into 40 GeV (top) and 130 GeV (bottom) binos. The pT of the hardest jet is 
plotted in the histograms (1 fb~^ luminosity). Matching is very important in the degenerate case 
when the contribution from initial state radiation is critical. The different colors indicate the 
contributions from Oj (orange), Ij (blue), and 2j (cyan). 
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To emphasize the importance of properly accounting for initial-state radiation using 
matching, Fig. [3] compares the Pt distribution for the hardest jet in a matched (left) and 
unmatched (right) dijet sample for a 150 GeV gluino directly decaying to a bino. The colors 
indicate the contributions from the different multiplicity samples: Oj (orange), Ij (blue), and 
2j (cyan). When the gluino-bino mass splitting is large enough to produce hard jets (top 
row), the Oj'' multiplicity sample is the main contributor. ISR is not important in this case 
and there is little difference between the matched and unmatched plots. The bottom row 
shows the results for a 130 GeV bino that is nearly degenerate with the gluino. In this case, 
only soft jets are produced in the decay and hard ISR jets are critical for having events pass 
the dijet cuts. Indeed, we see the dominance of the 2j multiplicity sample in the histogram 
of matched events. When ISR is important, the unmatched sample is clearly inadequate, 
with nearly 60% fewer events than the matched sample. 



B. Backgrounds 

The dominant backgrounds for jets + searches are /Z^ + jets, tt, and QCD. 
Additional background contributions come from single top and di-boson production (WW, 
WZ, ZZ), but these contributions are sub-dominant, so we do not consider them here. The 
missing transverse energy comes from vv and ^ l^v^ where the boson is 

produced directly or from the top quark. To reduce the background, a veto was placed 
on isolated leptons with > 10 GeV. However, these cuts do not completely eliminate the 
background because it is possible to miss either the electron or muon (or misidentify 
them). It should be noted that muon isolation cuts were not placed by PGS, but were applied 
by our analysis software. If the muon failed the isolation cut, then it was removed from the 
record and its four-momentum was added to that of the nearest jet. Additionally, the 
can decay into a hadronic r, which is identified as a jet. Because the DQ^ analysis did not 
veto on hadronic taus, we have treated all taus as jets in this study. 

QCD backgrounds can provide a significant source of low missing energy events, but are 
challenging to simulate. The backgrounds can arise from jet energy mismeasurement due 
to poorly instrumented regions of the detector (i.e., dead/hot calorimeter cells, jet punch- 
through, etc.). Additionally, there are many theoretical uncertainties - for example, in the 
PDFs, matrix elements, renormalisation, and factorisation/matching scales - that factor in 
the Monte Carlo simulations of the backgrounds. For heavy- flavor jets, there is the additional 

contribution coming from leptonic decays of the b-quarks. It is possible, for instance, to 
have the b-quark decay into a lepton and a neutrino, with the neutrino taking away a good 
portion of the b-quark's energy. Simulation of the QCD background is beyond the scope of 
Pythia and PGS and was not attempted in this work. To account for the QCD background, 
we imposed a tight lower bound on the of 100 GeV. Jet energy mismeasurement was 
accounted for by placing a lower bound of 90° and 50° on the azimuthal angle between the 
and the first and second hardest jets, respectively. In addition, an acoplanarity cut of 165° 
was placed between the two hardest jets. For the dijet case, the azimuthal angle between 
the and any jet with > 15 GeV and \vi\ < 2.5 was bounded from below by 40°. This 
cut was not placed on the threejet or multijet searches because of the greater jet multiplicity 
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in these cases. 

The W^/Z^ + nj and tt backgrounds were generated using MadGraph/MadEvent, with 
showering and hadronization in PYTHIA. PGS was again used as the detector simulator for 
jet clustering. The W^/Z^ backgrounds were matched up to 3 jets using the MLM match- 
ing procedure discussed in the previous section, with matching parameters = 10 GeV 
and Qmin = 15 GeV. The tt backgrounds were matched up to 2 jets with parameters 
Qmin = 14 GeV and = 20 GeV. For each of the separate backgrounds, 500K events 

were generated. The results approximately reproduce the shape and scale of the and 
Ht distributions published by the D(^ collaboration for 1 fb~^ [3j. In the dijet case, our 
results correspond to those of within ±20%. The correspondence is similar for the tt 
backgrounds in the threejet and multijet cases. For the W^/Z^ backgrounds, the correspon- 
dence is within ±30 — 40%. It is possible that this discrepancy is due to difficulties to fully 
populate the tails of the and Ht distributions with good statistics. In the case of the 
background, the modeling of the lepton detection efficiency in PGS might also play a 
role. Heavy ffavor jet contributions were found to contribute 2% to the W^/Z^ backgrounds, 
which is well below the uncertainties that arise from not having NLO calculations for these 
processes and from using PGS. 
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TABLE III: Differential cross section (in fb) for the monojet, dijet, threejet, and multijet samples 
of a theoretical model spectrum with a 340 GeV gluino decaying directly into a 100 GeV bino 
(4 fb~^). Some boxes show significant deviation from the signal limits shown in Table II: green 
indicates 0.5 < Xi ^ 2, blue indicates 2 < Xi < 3, and red indicates Xi > 3. All boxes with Xi > 1/2 
are included in the calculation of the total x^ value. 
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V. GLUINO EXCLUSION LIMITS 



A. No Cascade Decays 



For the remainder of the paper, we wiU discuss how model-independent jets + searches 
can be used to set hmits on the parameters in a particular theory. We will focus specifically 
on the case of pair-produced gluinos at the Tevatron and begin by considering the simplified 
scenario of a direct decay to the bino. The expected number of jets depends on the relative 
mass diflFerence between the gluino and bino. When the mass diflFerence is small, the decay 
jets are very soft and initial-state radiation is important; in this limit, the monojet search 
is best. When the mass diflFerence is large, the decay jets are hard and well-defined, so 
the multijet search is most effective. The dijet and threejet searches are important in the 
transition between these two limits. 

As an example, let us consider the model spectrum with a 340 GeV gluino decaying 
directly into a 100 GeV bino. In this case, the gluino is heavy and its mass diflFerence with 



the bino is relatively large, so we expect the multijet search to be most effective. Table III 



shows the differential cross section grids for the 1-4+ jet searches for this simulated signal 
point. The colors indicate the significance of the signal over the limits presented in Table II; 
the multijet search has the strongest excesses. 

Previously ^J, we obtained exclusion limits by optimizing the and Ht cuts, which 
involves simulating each mass point beforehand to determine which cuts are most appropri- 
ate. This is effectively like dealing with a 1 x 1 grid, for which a 95% exclusion corresponds 




500 



Gluino Mass (GeV) 



FIG. 4: The 95% exclusion region for at 4 fb"-*^ assuming 50% systematic error on background. 
The exclusion region for a directly decaying gluino is shown in light blue; the worst case scenario 
for the cascade decay is shown in dark blue. The dashed line represents the CMSSM points and 
the "X" is the current exclusion limit at 2 fb~^. 
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FIG. 5: 95% exclusion region (purple) for a 240 GeV gluino decaying into a bino through a wino. 
The dashed line is = + (9(m^o). The black dot at (m^, m^) = (60, 160), is the minimum 
bino mass for which a 240 GeV gluino is excluded for all wino masses. The inset shows the one-step 
cascade considered in the paper. 

to = 4. The approach considered here considers the significance of all such cuts, and only 
requires that a single n x n diflFerential cross section grid be produced for each search. 

Fig. [4] shows the 95% exclusion limit for directly decaying gluinos at 4 fb~^ luminosity 
and 50% systematic uncertainty on the background. The results show that such gluinos are 
completely excluded for masses below ~ 130 GeV. 

B. Cascade Decays 

In this section, we will discuss the exclusion limits for the decay chain illustrated in the 
inset of Fig. [5| In general, cascade decays are more challenging to see because they convert 
missing energy to visible energy.^ The number of jets per event is greater for cascading 
gluinos than directly decaying ones and the spectrum of jet energies depends on the ratio of 
gaugino masses. When rrig ^ rn^^^ two hard jets are produced in the decay of the wino to the 
bino. In the opposite limit, when ^ m^, two hard jets are produced in the decay of the 
gluino to the wino. When rrig < rriy^ < m^, four fairly hard jets are produced, diminishing 
the and making this region of parameter space the most challenging to see. In particular, 
the most difficult region to detect is when 

= ^B + C)(m^o). (12) 



^ For additional discussion of model-independent searches of cascade decays at the Tevatron, see [40] . 
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In the region of parameter space, where ^ m^^ the jets from the wino to bino decay 
become harder as the gauge bosons go on-sheU. 

Fig. [5] shows the values of and that are excluded up to 95% confidence for a 
240 GeV gluino (shaded region). The dark black dot, which represents the minimum bino 
mass for which a 240 GeV gluino is excluded for all wino masses, falls close to Eq. [12] (the 
dotted red line). 

The exclusion region in Fig. [5] is not symmetric about the line = + O(m^o). 
The asymmetry is a result of the hard lepton cuts. When the gluino and wino are nearly 
degenerate, the leptons from the gauge boson decays are energetic, and these events are 
eliminated by the tight lepton cuts, reducing the significance below the confidence limit. 
In the opposite limit, when the wino and bino are nearly degenerate, much less energy is 
transfered to the leptons and fewer signal events are cut. Additionally, the jets produced in 
this case are color octets and give rise to a greater number of soft jets, as compared to the 
singlet jets emitted in the gauge boson decays. The presence of many soft jets may decrease 
the lepton detection efficiency; as a result, it may be that even fewer events than expected 
are being cut. 

Figure [4] compares the 95% exclusion region for the cascade decay with that for the direct 
decay case. The "worst-possible" cascade scenario is plotted; that is, it is the maximum 
bino mass for which all wino masses are excluded. For the one-step cascade considered here, 
gluinos are completely excluded up to masses of ^ 125 GeV. 

C. t-channel squarks 

Thus far, it has been assumed that the squarks are heavy enough that they do not affect 
the production cross section of gluinos. If the squarks are not completely decoupled, they can 
contribute to t-channel diagrams in gluino pair-production. Figure [6] shows the production 
cross section for a 120 GeV (red), 240 GeV (blue), and 360 GeV (green) gluino, as a function 
of squark mass. When only one squark is light (and all the others are ^4 — 5 TeV), the 
production cross section is unaffected. However, when the squark masses are brought down 
close to the gluino mass, the production cross section decreases by as much as ^ 25%, 60%, 
and 75% for 120, 240, and 360 GeV gluinos, respectively. A reduction in the production cross 
section alters the exclusion region in the gluino-bino mass plane; while the overall shape of 
the exclusion region remains the same, its size scales with the production cross section [28]. 

It is worthwhile to note, however, that while the inclusion of squarks reduces the exclusion 
region for pair-produced gluinos by decreasing the production cross section, it also provides 
alternate discovery channels through gq or qq production. For example, if a gluino and 
squark are produced, with the gluino nearly degenerate with the bino, the subsequent decay 
of the squark will produce more visible energy than the gluino decay, thereby making the 
event more visible. 
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FIG. 6: Gluino production cross section as a function of squark mass: (red) rUg = 120 GeV, (blue) 
rrig = 240 GeV, and (green) rrig = 360 GeV. 

D. Monophoton Search 

Initial-state QCD radiation is important for gaining sensitivity to degenerate gluinos. 
Here, we will consider whether initial-state photon radiation may also be useful in the de- 
generate limit. Such events are characterized by small and a hard photon. 

The main benefit of the monophoton search is that the Standard Model backgrounds 
are better understood; unlike the monojet case, QCD is no longer an important background. 
Instead, the primary backgrounds come from processes such as Z^{^ vv) + 7, which is 
irreducible, and —> e^u where the electron is mistaken as a photon or W^{^ l^v) + 7, 
where the lepton is not detected. Other backgrounds may come from /Z^ + jet, where 
the jet is misidentified as a photon, or situations where muons or cosmic rays produce hard 
photons in the detector. 

The DQ^ Collaboration recently published results for their monophoton study, which 
searched for a Kaluza-Klein graviton produced along with a photon [41j. To reduce the 
Standard Model background, they required all events to have one photon with pt > 90 GeV 
and > 70 GeV. Events with muons or jets with pt > 15 GeV were rejected. They 
estimate the total number of background events to be 22.4 ± 2.5. 

To investigate the sensitivity of monophoton searches to degenerate spectra, we consider 
several points and compared them against 's background measurements. We considered 
several benchmark values for gluino and bino masses and did a simple cuts-based comparison 
between the monophoton search and an optimized monojet search. For example. Figure [5] 
shows that the monojet search safely excludes the case of a 140 GeV gluino and 130 GeV 
bino. A monophoton search (with the cuts used in the analysis) gives S/B = 0.48 and 
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S/^/B = 2.3 for this mass point; thus the monophoton search is roughly as sensitive but has 
a lower S/B value. Similarly, a 120 GeV gluino and 100 GeV bino is safely excluded by the 
monojet search, but the monophoton search only gives S/B = 0.39 and S/^/B = 1.86. 

There are several reasons why the monophoton search is not as successful as the monojet 
one. In the degenerate gluino region, the possibility of getting jets with a pt above the 
15 GeV threshold is significant (even though the mass difference is O{10 GeV)) because the 
gluinos are boosted. The monophoton search vetoes many events with such boosted decay 
jets. In addition, getting photon ISR is much more difficult than getting QCD ISR for 
several reasons - most importantly, because q^em <^ and because one is insensitive to the 
gluon-induced processes that contribute to the cross section. Despite these challenges, the 
significance of the monophoton search could still increase sensitivity. The monophoton does 
not fare significantly more poorly than the monojet one with the current set of cuts. Thus, 
it is possible that a more optimal set of cuts may increase the effectiveness of the search, 
especially given that the backgrounds are better understood in this case. Finally, the above 
estimates do not account for the photon detection rate in PGS, which may be different from 
that used by 's full detector simulator, from which the background estimates were taken. 



E. Leptons 

In this section, we address whether leptons from cascades can be used to augment the 
sensitivity of jets + searches. In the gluino cascade decay considered in this paper, it 
is possible to get leptons from the and boson decays. The 10 GeV lepton veto, 
however, eliminates most of these events. The exclusion limit for the gluino decay discussed 



in Sec. |VB| is not improved by removing the lepton veto; most of the irreducible backgrounds 
{W^ + nj and tt + nj) have a lepton and dominate over the signal when the veto is removed. 
The exclusion limit is not improved even if we require all events to have a certain number 
of leptons, or place cuts on lepton pt- 

The question still remains as to whether there is any region in parameter space where 
the jets + study places no exclusion, but a jets + + lepton study does. The lepton 
signal is useful for light gluinos (< 250 GeV) that are nearly degenerate with the wino. The 
signal point, a 210 GeV gluino decaying to a 50 GeV bino through a 170 GeV wino, is not 
excluded by the ordinary jets + analysis. We find here, though, that it has a significance"^ 
of 4: A for a pr cut of 50 GeV, but with a S/B 0.15. 

For high-mass gluinos, inclusion of the lepton signal does not increase the sensitivity 
of the search because the smaller production cross section decreases the signal significance. 
It might however be possible that lepton signatures are effective for high-mass gluinos in 
lepton-rich cascades that contain sleptons. Overall, though, these results indicate that while 
jets + + lepton searches may be useful in certain regions of parameter space, they should 
be combined with jets + searches to provide optimal coverage. 



^ Here, the estimate of the significance only accounts for the statistical error; it does not include the 
systematic uncertainty. 
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VI. CONCLUSION 



In this article, we discuss how model-independent bounds can be placed on the mass of 
the lightest color octet particle that is pair-produced at the Tevatron. The main aspects of 
the analysis focus on the advantage of running exclusive Ij — A^j searches, and placing limits 
using the measured differential cross section as a function of the visible and missing energy. 
We show that the exclusion reach can be significantly extended beyond those published by 
because the and Ht cuts used in their analysis were only optimized for points in 
CMSSM parameter space. The proposed analysis we present here opens up the searches 
to all regions of parameter-space, allowing us to set limits on all kinematically- accessible 
gluinos. We also show how the exclusion reach is degraded when gluino cascade decays are 
included, focusing on the example of an intermediate wino, which decays to the dark matter 
candidate. 

We have so far only focused on jet classification, and Ht as available handles for 
increasing the reach of jets + searches. However, in certain special cases, other techniques 
might be useful. For example, if the gluino decays dominantly to b jets, heavy fiavor tagging 
can be used advantageously. 

In our analysis of the cascade decays, we often found that the regions of highest signifi- 
cance in the diflFerential cross section plot were pressed down against the 100 GeV cutoflF in 
missing transverse energy. This lower limit was imposed to avoid regions where the QCD 
background dominates. If the 100 GeV limit could be reduced, then it would open up regions 
of high statistical significance that renders sensitivity to a larger region of parameter space. 
The numerous uncertainties in the theory and numerical generation of QCD events make it 
unlikely that precision QCD background will be generated in the near future. However, it 
may still be possible to reduce the cutoff by using event shape variables (i.e., sphericity). 

Looking forward to the LHC, jets + searches are still promising discovery channels 
for new physics. The general analysis presented in this paper can be taken forward to the 
LHC without any significant changes. The primary modification will be to optimize the 
jet Et used in the classification of the nj + searches. The backgrounds for the LHC 
are dominantly the same; however tt will be significantly larger and the size of the QCD 
background will also be different. Many of the existing proposals for searches at the LHC 
focus primarily on 4+^' + inclusive searches and are insensitive to compressed spectra; 
see [42j for further discussion on MSSM-specific compressed spectra at the LHC. By having 
exclusive searches over Ij + to 4+^ + the LHC will be sensitive to most beyond the 
Standard Model spectra that have viable dark matter candidates that appear in the decays 
of new strongly-produced particles, regardless of the spectrum. Additionally, having the 
differential cross section measurements will be useful in fitting models to any discoveries. 
Finally, it is necessary to confirm that there are no gaps in coverage between the LHC and 
Tevatron; in particular, if there is a light (~125 GeV) gluino, finding signal-poor control 
regions to measure the QCD background may be challenging. 
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